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We evaluate the effective potential for the conventional linear Walecka non perturbatively up to 
one loop. This quantity is then renormalized with a prescription which allows finite vacuum contri- 
butions to the three as well as four f PI Green's functions to survive. These terms, which are absent 
in the standard relativistic Hartree approximation, have a logarithmic energy scale dependence that 
can be tuned so as to mimic the effects of (f> 3 and (f> 4 type of terms present in the non linear Walecka 
model improving quantities such as the compressibility modulus and the effective nucleon mass, 
at saturation, by considering energy scales which are very close to the nucleon mass at vanishing 
density. 
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I. INTRODUCTION 

Quantum hadrodynamics (QHD) is an effective relativistic quantum field theory, based on mesons and baryons, 
which can be used at hadronic energy scales where the fundamental theory of strong interactions, quantum chro- 
modynamics (QCD), presents a highly nonlinear behavior. The Walecka model [l| to be considered here represents 
QHD by means of a lagrangian density formulated so as to describe nucleons interacting through the exchange of 
an isoscalar vector meson (u>) as well as of a scalar-isoscalar meson (0) which is introduced to simulate intermediate 
range attraction due to the s-wave isoscalar pion pairs. The original Walecka model (QHD-I) is described by the 
lagrangian density 



^ i C = ^ {id" - g v V") - (M - grf)] vb + \{d^d»4> - m 2 s 4> 2 ) - \Fu„F^ + \mlV^V" - U(cf>, V) + £ CT , (1.1) 
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l/~) ■ where ip, <p and ui denote respectively baryon, scalar and vector meson fields (with the latter being coupled to a 

conserved baryonic current). The term U(4>,V), which describes mesonic self interactions was set to zero in the 
iy~j . original model so as to minimize the many body effects while the term Cqt represents the counterterms needed to 
f^ ' eliminate any potential ultra violet divergences arising from vacuum computations. 

It is important to recall that, roughly, counterterms are composed by two distinct parts the first being a divergent 
piece which exactly eliminates the divergence resulting from the evaluation of a Green function at a given order in 
perturbation theory. The second piece is composed by a finite part which is arbitrary and can be fixed by choosing 
an appropriated renormalization scheme [2J. 

The important parameters are the ratios of coupling to masses, C 2 and C 2 , with Cf = {giM/rrii) 2 which are tuned 
to fit the saturation density p — 0.193fm~ 3 and binding energy per nucleon, BE = — 15.75MeV [l]. However, 
the QHD-I predictions for some other relevant static properties of nuclear matter do not agree well with the values 
quoted in the literature. For example, using the Mean Field Approximation (MFA) which considers only in medium 
contributions at the one loop level one obtains that, at saturation, the effective nucleon mass is M* &t ~ 0.56 M, 
which is somewhat low, while the compression modulus, K ~ 540 MeV, is too high. In principle, since this is a 
renormalizablc quantum field theory, vacuum contributions (and potential ultra violet divergences) can be properly 
treated yielding meaningful finite results. These contributions were first considered by Chin [3j|, at the one loop level, 
in the so called Relativistic Hartree Approximation (RHA) which produced a more reasonable value for the effective 
mass, M* at ~ 0.72 M. However, the compression modulus remained at a high value, K ~ 470 MeV. 

One could then try to improve the situation by also considering exchange contributions since both, MFA and RHA, 
consider only direct terms in a nonperturbative way. When vacuum contributions are neglected this approximation 
is known as the Hartree- Fock (HF) approximation producing M* at ~ 0.53 M and K ~ 585 MeV [l(. By comparing 
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the results from MFA, RHA and HF one sees how vacuum effects can improve the values of M* at and K. Then, the 
natural question is if the situation could be further improved by considering the vacuum in HF type of evaluations. 
The main concern now being the difficulty to deal with overall, nested, and overlapping type of divergences which 
surely arise due to the self consistent procedure. 

The complete evaluation of vacuum contributions within direct and exchange terms was performed by Furnstahl, 
Perry and Serot [4j (see Ref. [5j for an early attempt in which only the propagators have been fully renormalized ). 
This cumbersome calculation considers the nonperturbative evaluation and renormalization of the energy density up 
to the two loop level showing the non convergence of the loop expansion. Later, the situation has been addressed 
with the alternative Optimized Perturbation Theory (OPT) which allows for an easier manipulation of divergences 
[6|. Two loop contributions have been evaluated and renormalized in a perturbative fashion with nonperturbative 
results further generated via a variational criterion. However, saturation of nuclear matter could not be achieved with 
the results supporting those of Ref. [4j . 

Meanwhile, the compressibility modulus problem has been circumvented by introducing some more parameters, in 
the form of new couplings [7], to the original Walecka model. One then considers U(<j>, V) appearing in Eq. ()1.1[) as 

UW = ^ 3 + ^\ (1.2) 

where n — 2bMg 3 and A = 6cg 4 . This version of QHD is known as the nonlinear Walecka model (NLWM) and 
the main role of the two additional parameters, b and c, is to bring the compression modulus of nuclear matter and 
the nucleon effective mass under control. However, one may object to this course of action since the mesonic self 
interactions will increase the many body effects apart from increasing the parameter space. Notice also that terms 
like <j) n (n > 5) are not allowed since then, in 3+1 dimensions, one would need to introduce coupling parameters with 
negative mass dimensions spoiling the renormalizability of the original model, apart from increasing the parameter 
space. 

Heide and Rudaz [8j] have then realized that it is still possible to keep !/(</>, V) — while improving both K and 
M* at . The key ingredient in their approach is related to the complete evaluation (regularization and renormalization) 
of divergent vacuum contributions. Regularization is a formal way to isolate the divergences associated with a physical 
quantity for which many different prescriptions exist, e.g., sharp cut-off, Pauli-Villars, and Dimensional Regularization 
(DR). Within DR, which was used by Chin, one basically performs the evaluations in d— 2e dimensions taking e40 
at the end so that the ultra violet divergences show up as powers of 1/e. However, to keep the dimensionality right 
when doing d — > d — 2e one has to introduce arbitrary scales with dimensions of energy (A, or the related 1 A^jg). 
Chin has chosen a renormalization prescription in which the final results do not depend on the arbitrary energy scale 
while Heide and Rudaz chose one in which such a dependence remains, as in most QCD applications. Since the latter 
authors also worked at the one loop level their approximation became known as the Modified Relativistic Hartree 
Approximation (MRHA) and their main result was to show that it is possible to substantially improve K and M s * at 
by suitably fixing the energy scale, A. Moreover by choosing A = M the MRHA recovers RHA. In connection with 
neutron stars, the MRHA has been applied to the Walecka model in Refs. [9j. 

Here one of our goals is to treat the Walecka model using a formalism which is closely related to the one used in 
QCD and other modern quantum field theories. Within the QHD model, the temperature and density are usually 
introduced using the real time formalism employed in the original work of Walecka. Instead, we use Matsubara's 
imaginary time formalism treating the divergent integrals with DR adapted to the modified minimal subtraction 
renormalization scheme MS [2] which constitute the framework most commonly used within QCD. To obtain the 
ground state energy density, £, we will first evaluate the effective potential (or Landau's free energy), F, whose 
minimum gives the pressure, P. By choosing appropriate renormalization conditions we generate effective three- and 
four-body couplings, in J 7 , which are not present at the classical level. As we shall sec the numerical values of these 
effective couplings run with the energy scale, Aj^g, allowing for a good tuning of K and M* at which have their values 
improved at energy scales of about 0.92 M-0.98 M (M = 939 MeV) while the usual RHA results are retrieved for the 
choice A^jg = M. 

The MRHA proposed by Heide and Rudaz suggests that if one seeks to minimize many-body effects in nuclear matter 
at saturation, the choice A ~ M* at is the necessary one. Our philosophy is slightly different and perhaps simpler to 
implement. Since possible modifications in the behavior of K and M* seem to be dictated by the presence of K c ff0 3 
and X c s<p 4 type of terms we shall use the Chin- Walecka renormalization prescription to deal with <j> n (n = 0, 1,2) 
vacuum contributions by requiring that their respective contributions vanish at zero density (a requirement which 
was also adopted within the MRHA). However, as far as the vacuum contributions related to (j> 3 and 4> 4 are concerned 
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we advocate that one only needs to keep the finite energy scale dependent parts in the effective three- and four-body 
couplings. In this way, not only n e g and A e ff run with Aj^g but, as we shall see, one also retrieves the RHA at 
A M g = M._ 

Considering the effective potential at zero density we will choose an appropriate renormalization prescription for 
this particular model. Since our main goal is to improve K and M* at by quantically renormalizing k = — > K e s(A^g) 
and A = — > A cr f ( A^jg) we can keep M, m s , m v as representing the vacuum physical masses for simplicity. This choice 
means that, at kp = 0, all mass parameters (M,m s ,m v ) represent the effective vacuum masses and shall not run with 
Aj^g- as opposed to K er r and A c ff . In theories such as QCD the running of the couplings is dictated by the /3 function 
whose most important contributions come from the so-called leading logs , e.g. ln(Ajjg/M), which naturally arise in 
DR evaluations. The application of RG equations to the effective model Walecka model is beyond the scope of our 
work 2 . Nevertheless, our renormalization prescription to obtain a scale dependence so as to better control K and 
M* at is inspired by the leading logs role in the /3 function and the the renormalization scheme presented here proposes 
that one should preserve only the scale dependent leading logs which appear in the expressions for k c s4> 3 and A e ff</> 4 . 
As a byproduct, and contrary to the MRHA case, both quantities will display the same scale dependence. Here, this 
approximation will be called the Logarithmic Hartree Approximation (LHA). The numerical results show that the 
best LHA predictions for K and M* at are obtained at energy scales which are only about 5% smaller than that of the 
RHA, that is A^jg ~ 0.95 M. This is a nice feature since the values of the energy scale and that of the highest mass 
in the spectrum are almost the same whereas in the MRHA the optimum scale, set to be close to M* at is about 35% 
smaller than M. From the quantitative point of view, the LHA produces better results than the MRHA as will be 
shown. 

The work is presented as follows. In the next section the one loop free energy is evaluated using Matsubara's 
formalism. The renormalization of the vacuum contributions is discussed in Section III and the complete renormalized 
energy density is presented in Section IV. Numerical results and discussions appear in Section V where while our 
conclusions are presented in Section VI. For completeness, in the appendix, we discuss a case in which m s does not 
represent the physical mass. 

II. THE FREE ENERGY TO ONE LOOP 

In quantum field theories the effective potential (or Landau's free energy) , _F, is defined as the generator of all one 
particle irreducible (1PI) Green's functions with zero external momentum. The standard textbook definition (for one 
field, 4>) reads 

OO 

^c) = £f (n) (0)#, (2.1) 

n=0 

where have absorbed non relevant factors of I and n\ by defining r(")(0) = (— i) n r(")(0)/n! with r(™)(0) representing 
the 1PI n-point Green's function and <j) c representing the classical (space-time independent) scalar field. In practice, 
this quantity incorporates quantum (or radiative) corrections to the classical potential which appears in the original 
lagrangian density. While the latter is always finite the former can diverge due to the evaluation of momentum integrals 
present in the Feynman loops. One way to obtain this free energy density is to perform a functional integration over 
the fermionic fields [2J. To one loop this leads to 

H<f>o V c ) = -^fv^Vf + ^4>l + ij -0^trln[r(k, - 9v V c>li ) (M g s <t> c )] . (2.2) 

Notice that this free energy density contains the classical potential (zero loop or tree level term) present in the 
lagrangian density plus a one loop quantum (radiative) correction represented by the third term. Working in the rest 
frame of nuclear matter we assume that the classical fields are time-like (V c ^ — 8^^V C ^). Then, after taking the trace 
one can write the free energy as 

H4>c V c , ) = -^K 2 + ^ + i 7 1 (0j ln[-(*o - g v V cfi f + k 2 + (M 5s c ) 2 ] , (2.3) 



See Ref. |10ll for a RG investigation of the Walecka model. 



where 7 = 4(2) is the spin-isospin degeneracy for nuclear (neutron) matter. To obtain finite density results one may 
use Matsubara's imaginary time formalism with k Q — > i(u> n ~ *m) where /i represents the chemical potential while, for 
fermions, us n — (2n + \)itT (n = 0, 1, ...) is the Matsubara frequency with T representing the temperature. Then, the 
free energy reads 



^c,Vc,0) = -^ + ^ 



^ T E J j^f ln{K " ^ ~ 9vVcfi)]2 + k2 + (M ~ 9s<j)c)2] 
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The Matsubara's sums can be performed using 
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where E 2 (k) — k 2 + (M — g s (j> c ) 2 and /jf = /j, — g v V C M- Being interested in the T — case one may take the zero 
temperature limit of Eq. (|2.5|) which is given by 
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Then, at T = and \i ^ 0, the one loop free energy for the Walecka model becomes 

™2 , d 3 k ^ 
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Power counting shows that A(</> c ) is a divergent quantity while the [i dependent term of Eq. (|2.7[) is convergent due 
to the Heaviside step function. 



III. THE RENORMALIZED VACUUM CORRECTION TERM 



In order to renormalize the vacuum correction term one must first isolate the divergences which is formally achieved 
by regularizing the divergent integral. Here we use DR performing the divergent integrals in 2u = 3 — 2e dimensions. 
Then, in order to introduce the MS energy scale, Aj^g, commonly used within QCD one redefines the integral measure 
as 



d 3 k 

(2tt) 3 



-> 



47T 



e/2 



d 2uJ k 

(2tt) 2 " 



(3.1) 



where je = —0.5772... represents the Euler-Mascheroni constant. Note that, with this definition, irrelevant factors 
of je and 4n are automatically cancelled but the results of Refs. [3, @, @| can be readily reproduced by using 
Aj^g = \/47re -TB A. The integral can then be performed yielding |2| 
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As one can see, by expanding the the term proportional to 1/e, there are five potentially divergent contributions 
ranging from g° to g 4 while all terms of order g n (n > 5) are convergent. The divergent terms proportional to T^ n '(j)™ 
(11 = 0, ...,4) are respectively 
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The counterterms contained in Cqt needed to render the free energy finite are |3j, | 
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where the a n coefficients have the general form 



a n ~ 5 S 



+ /n( A ls) 



(3.9) 



Now, within the MS renormalization scheme generally adopted within QCD one simply sets /„ = and the 
counterterms have only the bare bones needed to eliminate the 1/e poles while the final finite contributions depend 
on the arbitrary energy scale. If one adopts this scheme within the Walecka model the free energy would look like the 
dashed curve in Fig. Q] which shows T versus C for the values 3 A^jg = 0.9 GeV, M — 1 GeV, m s — 0.55 GeV, and 
g s = 1. As it is well known, within this scheme M, m s , and m v do not represent the measurable physical vacuum 
masses which are instead taken as mass parameters whose values, like the values of the couplings, run with A^g in a 
way ultimateley dictated by RG equation. 
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FIG. 1: (color online) The free energy, in the c/> c direction, as a function of the classical field for A^jg = 0.9 GeV, 

M = 1 GeV, m s = 0.55 GeV, and g s = 1. The dashed line is the MS renormalization scheme result. The dotted- 

dashed corresponds to T^ ) = rW = while MS is used in the remaining three 1PI function. The same situation but 

with r^ 2 ) = is represented by the dotted line. The continuous line represents the LHA prescription. 



3 Note that some of the these values are close to the ones which will later be used in our numerical procedure. However, at this stage 
they are not intended to represent any realistic physical situation apart from letting us compare possible different shapes of T. 



If instead, like Chin, one adopts the so-called on-mass renormalization scheme the counterterms completely eliminate 
the total contributions represented by Eqs (j3.3H3.7p . Within this choice the results are scale independent while M, m s , 
and m v represent the measurable physical masses at zero density whereas the three and four-body mesonic couplings 
vanish in agreement with the tree level result displayed by the original lagrangian density. The free energy generated 
by this scheme is represented by the dashed line in Fig. [2j Considering the relevant Kp = case, let us find a hybrid 
alternative scheme between the MS and the on-mass-shell so that a residual, scale dependent, contribution survives 
within the three and four 1PI Green's function given by Eqs (|3.6p and (|3.7[) . 
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FIG. 2: (color online) The free energy, in the </> c direction, as a function of the classical field for A^jg = 0.9 GeV, 
M = 1 GeV, m s = 0.55 GeV, and g s = 1. All curves represent the LHA prescription at the different scales A^jg = 
0.9 GeV < M (continuous line), A^g = 1 GeV = M (dashed line) and A^jg = 1.1 GeV > M. The dashed line also 

corresponds to the RHA result. 



To do that, let us analyze each of the arbitrary /„ terms contained in the counterterm coefficients from the physical 
point of view starting with /q which is contained in the field independent r^ '. This contribution is renormalized by 
the constant counterterm «o which can be referred to as the "cosmological constant" [ll| . In practice, the only effect 
this term has is to give the zero point energy value and by its complete elimination one assures that J-(4> c = 0) = 
which, within the Walecka model, will later assure that the pressure as well as the energy density vanish at kp = 0. 
Therefore, as in the on-mass shell prescription, we can impose that /o be exactly equal to the finite part of the T^ 
term. It is important to point out that even if one uses the MS scheme this term can be absorbed in a vacuum 
expectation value subtraction of the zero point energy so that the exact way in which it done is not too relevant for 
the present purposes. 

The effect of the the linear (tadpole) term T^ l '<f> c is to shift the origin so that the minimum is not at the origin 
(4> c ^ 0) as shown by the dashed line of Fig. Q] Also any finite contribution left in the tadpole will cause direct 
terms to contribute to the baryon self energy which, at the present level of approximation, means that M does not 
represent the physical nucleon mass at kp = 0, M* ac . This can be understood by recalling that the baryon self-energy 
is Tib ~ (gs/mjJF^fAijg) so that the vacuum effective baryon mass is given by M* ac = M + Ss(Aj^g) and since 
M* at = 939 MeV one sees that M, as well as g s and m s , should depend on Aj^g. However, for the purposes of 
controlling K and Af s * at the renormalization of the baryonic vacuum mass from M to M* ac does not generate the 
wanted 4> 3 and 4> 4 vertices. Therefore, for simplicity, we can also set /i so as to completely eliminate the tadpole 
vacuum contribution. This choice for /o and /i together with f 2 = f% = f& = produces the dot-dashed line of figure 
[TJ The term T^ 2 ) represents a (momentum independent) vacuum correction to the scalar meson mass, m s . As in the 
previous case, getting rid of this term assures that m s be taken as the physical mass simplifying the calculations since 
( TO s vac) 2 = m s + ^s ( Ajyig) where E s (Aj^g) ~ f( 2 )(Aj^g). Fixing f 2 so as to completely eliminate the r^ 2 ) contribution 
produces the dotted line of figure [1] In summary, so far we have adopted the usual Chin- Walecka on-mass shell 
renormalization conditions for /i, /2, and fg, so that: the vacuum energy is normalized to zero, 4> c = is the minimum 
of J- (also meaning that M = M* ac ), while m s represents the vacuum scalar meson mass. In this approach, none 
of the vacuum mass parameters present in the original lagrangian density run with Aj^g. Note that, physically, our 
choice was also inspired by the NLWM observation that the compressibility modulus is improved by the introduction 
of 4> 3 and (j) 4 terms which is consistent with our choice of neglecting any corrections to terms proportional to (/n/ji/; 
and (j) 2 which are directly related with the scalar meson and baryon masses. 

Now, taking fy=0 and fi = would leave us with the wanted </> 3 and </) 4 scale dependent terms. However, inspection 
of Eq. p.6p and Eq. (|3.7p shows that these contributions would vanish at different scales, given by Ajjg = Me 1 ' 3 



and Aj^g = M e 4 / 3 respectively. As already emphasized the NLWM controls the compression modulus with the kcj> 3 
and X(f> terms so one can impose that, within our approach, both K e & and A c g arise at the same energy scale. This 
can be achieved by imposing that r' 3 ) = I^ 4 ) = at A^jg = M in which case the RHA is always reproduced. Finally, 
in our RG-NLWM inspired prescription we also impose that any A^jg dependence should be left within the leading 
logs which naturally emerge within DR, as shown by Eqs (j3.3H3.7p , and which are the main contributing terms to the 
P function. In this case, the 0:3 and ay counterterms also eliminate the Aj^g- independent constants in Eqs (|3.6j) and 
(13. 7[) . One then obtains the continuous line of Fig. [TJ The complete finite, scale dependent, vacuum contribution is 
then given by 
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The free energy obtained with this finite vacuum contribution term is shown in Fig [2] for Aj^g- < M (continuous 
line), Ajyfg > M (dot-dashed line) as well as for A^jg = M (dashed line) in which case the usual RHA is retrieved. As 
one can check, the first term in Eq. (|3.10[) is just the RHA vacuum correction Q so that, in view of Eq. (|1.2j) . one 
can write 
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where k cS = 2g->Mb c s(A MS ) and A e ff = 6g 4 c e ff (A^g) with 
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and & c ff(Aj^jg) = — 3c ff(Aj^g). In this way not only n e s and A e ff vanish at the same scale but an inversion of their 
respective signs happen at the same time. We have then achieved our goal by quanticaly inducing k = — > K c g(A^jg) 
and A = — > A c ff(A^|g-) in a way that all the scale dependence is contained in the leading logs and also achieving 
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For comparison purposes let us quote the MRHA result 
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One notices that the major difference between the MRHA and our prescription amounts to the finite contribution 
contained within the cubic term where the scale dependence is not restricted to the leading log being also contained 
in an extra linear term which does not naturally arise when expands the DR results for the loop integrals in powers 
of e, as shown by Eqs. (j3.3H3.7p . 

IV. RENORMALIZED ENERGY DENSITY 

To obtain the thermodynamical potential, £7, one minimizes the free energy (or effective potential) with respect to 
the fields. That is, il — J-{<t Ci Vq) = —P, where P represents the pressure. Then, the LHA renormalized pressure is 
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where E* — (k 2 + M *) 1 / 2 with M* = M — g s cj) c while the Fermi momentum is given by k 2 F = (/J, — g v V 0>c ) 2 — M* 
For the vector field one gets 
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where ps = (jkp)/(6ir 2 ) is the baryonic density whereas for the scalar field the result is 
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represents the scalar density and 
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To get the energy density, £ , one can use the relation £ = — P + pps obtaining 
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M 



A 
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V. NUMERICAL RESULTS 

Let us now investigate the numerical results furnished by LHA for the baryon mass at saturation as well as for the 
compressibility modulus, with the latter given by 



K 



dk 2 \p B ) 



k—kf 



Pb 



d 2 



dp 2 B \pb) 



(5.1) 



PB=Pa 



Table U shows the coupling constants and saturation properties for some values of the renormalization scale (A^jg) 
that yield BE — —15.75 MeV and fcjr = 1.42 fm _1 (280.20 MeV). These values are chosen just in order to compare 



TABLE I: Coupling constants and saturation properties for some values of the renormalization scale (Ajjg-) that yield 
BE = -15.75 (MeV) and k F = 1.42 fin -1 . The meson masses are m s = 512 MeV and m v = 783 MeV. 



^ms/M K (MeV) M^/M 



C ,. 



C: 



<-ii 



Hqb/M 



(MS) 


1.030 


1279.408 


0.606 


171.339 


176.984 


119.138 


52.619 


-6.859 


49.753 


(MS) 


1.020 


910.234 


0.646 


151.744 


184.093 


105.512 


54.736 


-4.875 


36.063 


(MS) 


1.010 


639.833 


0.684 


132.232 


185.875 


91.945 


55.263 


-2.485 


18.474 


(MS) 


1.005 


542.279 


0.702 


123.626 


185.609 


85.961 


55.184 


-1.243 


9.233 


(MS) 


1.000 


468.140 


0.718 


114.740 


183.300 


79.782 


54.497 


0.000 


0.000 


(MS) 


0.990 


371.437 


0.745 


99.784 


177.933 


69.383 


52.901 


2.351 


-17.099 


(MS) 


0.980 


314.086 


0.767 


88.623 


173.525 


61.622 


51.591 


4.551 


-32.689 


(MS) 


0.975 


294.260 


0.776 


84.456 


172.456 


58.725 


51.273 


5.651 


-40.463 


(MS) 


0.970 


277.989 


0.784 


80.307 


170.683 


55.840 


50.746 


6.694 


-47.684 


(MS) 


0.960 


253.249 


0.798 


73.691 


168.648 


51.240 


50.141 


8.811 


-62.391 


(MS) 


0.950 


235.660 


0.809 


67.923 


166.440 


47.229 


49.484 


10.855 


-76.356 


(MS) 


0.940 


222.493 


0.818 


63.025 


164.576 


43.824 


48.930 


12.875 


-90.058 


(MS) 


0.920 


202.507 


0.832 


55.411 


162.702 


38.529 


48.373 


17.054 


-118.611 


(MS) 


0.900 


188.175 


0.843 


49.467 


161.826 


34.396 


48.112 


21.375 


-148.267 


(MS) 


0.8595 


166.351 


0.8595 


40.936 


164.038 


28.464 


48.770 


31.349 


-218.926 


(MS) 


0.850 


162.638 


0.863 


39.348 


165.080 


27.360 


49.080 


33.971 


-237.992 


(MS) 


0.800 


144.532 


0.876 


32.511 


172.680 


22.606 


51.340 


49.902 


-357.552 


(MS) 


0.700 


118.409 


0.893 


23.328 


200.551 


16.221 


59.626 


99.833 


-770.891 


(MS) 


0.600 


98.285 


0.905 


17.052 


256.570 


11.857 


76.281 


206.893 


-1806.98 


(MS) 


0.500 


81.019 


0.914 


12.239 


397.387 


8.510 


118.147 


541.142 


-5881.97 


(MS) 


0.400 


65.319 


0.921 


8.235 


1290.240 


5.726 


383.600 


4185.070 


-81967.5 


(RHA) 


_ 


468.140 


0.718 


114.740 


183.300 


79.782 


54.497 


_ 


_ 


(MFT) 


- 


546.610 


0.556 


195.900 


267.100 


136.210 


79.423 


- 


" 



with the original Walecka Model (QHD-I) [1]. The meson masses are m s = 512 MeV and m v — 783 MeV. This table 
shows that some of the best LHA values are obtained with A^jg- values which are very close to M. Since at Aj^g- = M 
the RHA result is reproduced one concludes, based on our results, that slight decrease from the RHA energy scale 
produces an enormous effect on the values of both, K and M s * at . 

Figures [3] (a) and (b) show the binding energy per baryon, BE = E/A — M, and the effective baryon mass as 
functions of the Fermi momentum for some A^jg- values, shown in table [I] One easily sees the effect of considering 
the vacuum contribution and its improvements on the compressibility and the effective mass. As expected, when 
A^jg = M, the RHA results are recovered. From figures |4] (a) and (b) it is possible to see some properties obtained 
in table U within the LHA approach, as functions of Aj^g/M. One notes from figure |4] (a) that when A^jg- increases 
the value of the nuclear compressibility (K) also increases and M* at decreases. The crossing point in figure @] (a) 
represents the RHA values of K and M* which occurs when we set A^g- = M. Figure [4] (b) shows the effective 
couplings that arise due to the LHA as functions of A^g/M. Similarly when Aj^g- reaches the value M the RHA 
results are recovered and the effective couplings vanish. 

To compare our numerical results with those provided by the MRHA let us make a remark concerning the effective 
nucleon mass. From a non-relativistic analysis of scattering of neutron-Pb nuclei it has been found [L2| that M* at /M ss 
0.74 to 0.82 which can be viewed as approximately describing the Landau effective mass [13fl . The relativistic isoscalar 
component known as the effective mass defined in Eq. (|4.8[) can be called the Dirac effective mass and is related to 
the Landau effective mass. Therefore, the range expected for the Dirac effective mass at saturation density lies in 
the range M* at /M sa 0.70 to 0.80 whereas for the nuclear compressibility at saturation the most widely accepted 
values are K w 200 MEV to 300 MeV [14]. For this range of K and according to table [III the MRHA predicts 
M* at /M « 0.80 to 0.85 for K/M « 1.185 to 1.466. However, one should note that this MRHA energy scale range is 
not unique and can also be reproduced with A/M sw 0.753 to 0.778 which in turn leads to a rather low range for M* at 
values, M* at /M ps 0.65 to 0.69. Our results, shown in tables HI and [Til seem to produce a better agreement for this 
range of K giving the unique range M* at /M sa 0.76 to 0.83 for A^g/M « 0.920 to 0.977 with k cS > and A off < 0. 
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FIG. 3: (a) Binding energy per nucleon as a function of the Fermi momentum for different values of our scale Aj^r. 
We also plot the MFT and RHA results for comparison purposes. The saturation properties are: BE = —15.75 MeV 
and kp — 1.42 fm -1 (280.20 MeV). (b) Similar as figure (a) but for the effective baryon mass M* as a function of 



the Fermi momentum for different values of the scale. Note that when A 



MS 



1 we reproduce the RHA results. 
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FIG. 4: Some properties obtained in table U within the LHA approach, as a function of the renormalization scale 
Aj^g in units of the baryon mass, (a) Compression modulus and effective baryon mass x Ajm/M and (b) effective 
couplings that arise due to the LHA as functions of A^jg/M. They vanish when A^-^/M — 1 and the RHA results 

are reproduced as one sees in table |U 



As a last remark we would like to point out that if one chooses A^a/M — 0.9805, (g s /m s ) 2 = 9.468 fm 2 and 



(W T 



4.879 fm 2 the LHA approach reproduces the same saturation properties as performed by the so-called 



GM2 parameter set according to [15j: K = 300 MeV, M* at /M = 0.78, BE = -16.3 MeV, k F = 1.313 fm" 1 and 
po = 0.153 fm -3 . The resulting effective couplings are: & e ff = 0.005986 and c e g = —0.001995. 

In the appendix we show that leaving a leading log dependence also in the two point Green's function, T^ 2 ', only 
increases the numerical complexity without producing results better than the ones generated by the simplest LHA 
version employed so far. 



VI. CONCLUSIONS 



We considered the simplest form of the Walecka model to analyze how the values of the compressibility modulus as 
well as the baryon mass, at saturation, can be improved by adopting an appropriate renormalization scheme in which 
cubic and quartic effective couplings are radiatively generated. With this aim we have evaluated the effective potential 
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TABLE II: Comparison between the LHA and the MRHA approaches with other estimates. 
A s /A/ K (MeV) M* /M 

MRHA [8] 1.185 - 1.466 300 - 200 0.80 - 0.85 

MRHA [8] 0.778 - 0.753 300 - 200 0.65 - 0.69 

LHA 0.977 - 0.920 300 - 200 0.76 - 0.83 

Estimates [12-14] - 300 - 200 0.70 - 0.80 

Where A s = A for MRHA and LHA is given by: A s = A^jg . 

to the one loop level using Matsubara's formalism to introduce the density dependence. The vacuum contributions 
have been evaluated using dimensional renormalization compatible with the MS renormalization scheme. 

We have then chosen the renormalization conditions in such a way so that all the mass parameters appearing in the 
original lagrangian density represent the physical mass at zero density and therefore do not run with the energy scale. 
For our purposes the most important part was to renormalize the values of the cubic and quartic terms (kc/) 3 /3! and 
A</> 4 /4!) which vanish at the classical (tree) level in the original model. We have then allowed only scale dependent 
logarithms, which naturally arise within DR, to be present in the final finite expressions and, contrary to the MRHA 
prescription, we obtained that both couplings have exactly the same type of scale dependence. In other words, the 
parameters b and c contained in the cubic and quartic terms have been dressed by one loop vacuum contributions so 
that b = ->• b cS = 3/tt 2 ln(Aj^/M) and c = ->• — & cff /3. 

In our approach each value of the energy scale produces only one value for K and M* at while two values can be 
obtained within the MRHA. In our case the best values for these physical quantities occur at energy scales very close 
to the highest mass value, M . Since the RHA is obtained for Aj^r = M one concludes that a small variation around 
this value of the energy scale can significantly improve both K and M* at as shown by our numerical results which 
predict M* at /M « 0.76 to 0.83 and K « 200 MeV to 300 MeV at A^g/M sw 0.920 to 0.977. These results turn to be 
in excellent agreement with the most quoted estimates M* at /M sw 0.70 to 0.80 and K sw 200 MeV to 300 MeV. To 
achieve these K values the MRHA predicts either M* a JM « 0.80 to 0.85 or M* at /M s=s 0.65 to 0.69 in the two possible 
energy scale ranges. Recalling that at A-^/M = 1 the (RHA) results are M* at /M = 0.718 and K = 468.14 MeV 
one may further appreciate how a small tuning of the energy scale within the LHA greatly improves the situation. 
To compare the LHA with the MRHA we recall that the philosophy within the latter is that many-body effects in 
nuclear matter at saturation can be minimized by choosing the energy scale close to M* at in which case the values 
M* at /M = 0.731 and K = 162 MeV are reproduced. Although the former seems reasonable the latter seems too low 
according to the above quoted estimates. The philosophy of the LHA, proposed in the present work, is to keep only 
the scale dependent leading logs in the finite parts of the effective cubic and quartic couplings. 

In practice, the main difference between the two approximations is reflected by the fact that the MRHA effective 
cubic coupling, apart from the logarithmic term, also displays a term which depends linearly on the energy scale 
accounting for the numerical differences cited above. It is worth pointing out that, within the LHA as well as the 
MRHA, a given scale sets both K c g and X c g so that both K and M* at cannot be separately tuned as in the NLWM 
where n and A can be set separately. However, even in an effective theory such as the Walecka model, an increase in 
the parameter space as the one generated by the NLWM can be viewed as an unwanted feature and the LHA succeeds 
in improving the values of K of M* at without the drawback of increasing many body effects and parameter space. 
The method proposed here should be easy to be implemented within many existing MFA or RHA applications where 
A LHA can be added to the energy density (in the MFA case) or used to replace the existing A RHA in a RHA type of 
calculation. 

In principle the LHA philosophy could be extended to the two loop level in a calculation similar to the one performed 
in Refs. [J] and [6j. Then by tuning the energy scale appropriately one could try to reduce the size of the two loop 
corrections producing physically meaningful results. 

Our method can be extended to applications related to neutron stars and evaluation of other physical quantities, 
such as the symmetry energy. In particular, models and/or approximations which lead to low effective masses at 
saturation are not suitable for neutron stars calculations since as the density increases the effective mass vanishes so 
quickly that higher densities cannot be properly reached as needed [16j| . In principle, the LHA has potential to correct 
this problem without the need to introduce extra mesonic interactions with their respective parameters. 
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Appendix 

For completeness, let us check numerically the effects of leaving a leading log dependence also in the two point 
Green's function with zero external momentum, T^ 2 ' , given by Eq. (J3.5I) . Then, 



,LHA; 



(M*, A^) -> A^ A (M*,A m ) - (M - M*YM 



2 j\ /rl I 



16tt 2 



6 In 



M 



(A.l) 



and 



A 



/LHA 



{M*,k m ) 



A /LHA 



(M*,A- 



sis<-(M-M*)M 2 -L 



6 In 



M 

A~~ 



MS 



(A.2) 



In this case, the effective potential gives a first (momentum independent) correction to the effective scalar mass in 
the vacuum is m* vac . Then, for each energy scale, apart from the BE requirement one also has to fix the parameter 
set so that the effective scalar meson mass, in the vacuum m* vac = 512 MeV. This effective mass is obtained by 
considering one loop momentum independent self energy 



I * \2 2 11^2 2 7 

(TO S ,vac = m s - M 9sl^2 



6 In 



M 



A 



MS 



(A.3) 



which clearly indicates that m s (as well as g s ) must run with the energy scale. However, this more cumbersome 
approach has almost no effect in our best results for K and M* at as table IIIII shows indicating the adequacy of the 
LHA simple prescription previously adopted. 



TABLE III: Same as in table Ufor the case in which r^ 2 ) has a non vanishing leading log and m s runs with the energy 

scale. 



A m /M K (MeV) M s * at /Af 



(-'v 



ci 



gi m s (MeV) 



(MS) 


1.000 


468.140 


0.718 


(MS) 


0.975 


294.260 


0.776 


(MS) 


0.950 


235.660 


0.809 


(MS) 


0.920 


202.507 


0.832 



114.740 183.300 79.782 

84.456 74.106 58.725 

67.923 46.297 47.229 

55.411 31.755 38.529 



54.497 512.000 

22.032 641.594 

13.765 671.839 

9.441 687.840 
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